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FOREWORD 

This Interim report presents the results of work performed under 
Contract NAS8-34975 for the National Aeronautics aid Space Administration, 
George C. Marshall Space Flight Center, Huntsville, Alabama. This work was 
performed by personnel In the Product Engineering & Development Section of 
the Lockheed-Huntsvllle Research & Engineering Center and by the Computa- 
tional Mechanics Company, Austin, Texas, subcontractor to Lockheed during 
this effort. 

The period of performance for this study was from August 1982 through 
August 1983. The MSFC Contracting Officer's Representative for this study 
is Larry A. Kief ling, ED22. 
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SUMMARY 

This report contains the results of a study whose objective was to add 
the capability of analyzing a coupled dynaalc system of flowing fluid and 
elastic structure to the SPAR computer code. A comprehensive literature re- 
view was performed and a method was developed and adopted for use In SPAR. 
The method utilizes the existing assiimed-stress hybrid Plan element cur- 
rently In SPAR. An operational module was Incorporated In SPAR which pro- 
vides the capability for analyzing the flow of a two-dimensional, Incompres- 
sible, viscous fluid within rigid boundaries. Equations were developed to 
provide for the eventual analysis of the Interaction of such fluids with an 
elastic solid. 
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1. INTRODUCTION 

Virtually every structure Is In contact with a fluid, be It air, water, 
or a gas flowing by design or by Its natural course over and through the 
structure's surfaces. The fluid thereby exerts loads on the structure pro- 
ducing deforma *:lons which may. In turn, alter the flow of the fluid. 

In most situations encountered in the design of engineering systems, 
this fluid-structure li>.:eractlon is Insignificant, and the structure and the 
fluid can be analyzed Independently. There are Important cases, however, in 
which the Interaction of fluid and structural behavior is an Intrinsic fea- 
ture of the response of both media, and this Interaction must be taken into 
account In any rational analysis and design. Such is the case, for example, 
in the analysis of flutter phenomena In aircraft, the sloshing of fuels or 
other liquids in flexible tanks, flow-induced vibrations of sulnaerged struc- 
tures or tall buildings, the safety analysis of nuclear reactor components - 
particularly the study of pressurized reactor cores - the flow of liquids In 
flexible pipes as in the flow of blood in elastic arteries or oil or water 
In rubber hoses, the effects of underwater explosions on submerged struc- 
tures, etc. 

Fluid-structure interaction problems such as these are inherently non- 
linear: the Homaln of fluid media obviously changes with the deformation of 

the structure and pressures exerted by the fluid act on material surfaces 
the locations of which depend upon the deformation. There are, however, 
significant classes of fluid-structure interaction problems for which useful 
results can be obtained by using only linearized equations. Indeed, the 
bulk of the work published on this subject deals with one special case or 
another for which the analyses can be dealt with using linear or mildly 
nonlinear theories. 
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In very recent times, important applications have arisen In which a 
study of rather general and highly nonlinear fluid-structure Interaction 
phenomena Is needed. Because of the formidable mathematical difficulties 
Inherent In such nonlinear problems, most analysis procedures In use today 
are designed for cominiter Implementation. Indeed, for mere than two dec-* 
ades, a slgnlflrsnt volume of literature on the numerical analysis of fluid- 
structure interaction problems has accumulated and much of the work over the 
last decade has Involved the development of finite element methods and has 
primarily focused on problems of nuclear reactor safety. 

This report contains a survey and a critical analysis of current nu- 
merical schemes used for fluid-structure Interaction problems, bpeclal 
emphasis Is placed on finite element methods and in various models and algo- 
rithms now In use or under study for a wide class of such problems. Ue will 
adopt a deductive approach to this subject, considering first the formula- 
tion of very general models for fluid-structure Interaction and then reduc- 
ing these to various special cases that may arise In specific applications. 

Following this Introduction, Section 2 contains a brief survey of some 
of the relevant literature. The principal sources consulted In the prepara- 
tion of this document are collected in the Bibliography. In Section 3, we 
discuss so-called mixed Eulerlan-T.agrangian descriptions of motion and the 
corresponding klnematlcal equations. An attempt Is made to present this 
subject in a relatively thorough and complete way and to provide some clar- 
ity and precision in deriving fundamental klnematlcal relations that are 
critical to subsequent developments. Derivation of the equations of motion 
of an arbitrary fluid or solid in such a mixed reference frame Is taken up 
In Section 4. These equations provide the basis for derivation of the seml- 
dlscrete systems governing finite-element models of fluid structure Interac- 
tion, discussed In Sections 5 and 6. Interface conditions are ttken up In 
Section 7, Section 8 contains a derivation of a linear fluid-structure 
Interaction model for nearly Incompressible viscous flows. Modifications to 
the SPAR code to permit Incompressible viscous flow calculations are given 
in Section 9. Results of example problems arc presented In Section 10. 
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2. FINITE ELEKENT METHODS FOR FLUID-STRUCTURE INTEiUCTION 
2.1 GENERAL REMARKS 

Interest in the development of finite element methods for fluid- 
structure Interaction problems began primarily In the mid-1970s when concern 
over the structural Integrity of nuclear containment vessels called for 
better methods for modeling the reactor core, the liquid coolant, surround- 
ing gases and the vessel walls under various conditions. The names 
Belytschko, Donee, Kennedy, and Liu are prominent In this body of litera- 
ture, and several surveys of literature on computational methods for fluid- 
structure Interaction In nuclear reactors have been presented by Belytschko 
and Donea and their associates. In this regard, see, for example, 

Belytschko (Refs. 1 and 2), Kennedy and Belytschko (Refs. 3 through 5), 
Belytschko and Kennedy (Refs. 6 and 7), and the references therein to the 
series of articles by Donea, Fasoll-Stella, and Giuliani (Refs 8 through 
10), Donea, Giuliani, and Halleux (Ref. 11), Donea (Ref. 12), and the dis- 
sertation of Liu (Ref. 13). The voluminous collections of proceedings of 
the biannual SMIRT (Structural Mechanics in Reactor Technology) conferences 
contain numerous papers on computational methods for fluid-structure Inter- 
action problems and there one can find a heavy emphasis on finite element 
methods. 

A significant but relatively smaller collection of papers ^las been pub- 
lished on finite element methods for flow-induced vibrations o: structures, 
wave effects on submerged structures, and sloshing of liquids In elastic 
tanks. He shall cite representative examples of this literature later. 
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2.2 VARIOUS DESCRIPTIONS OF MOTION 


The first end nost fundsnental issue that confronts one in modeling 
fluid-structure Interaction is the choice of an appropriate framework for 
the description of the motion. Traditionally, in solid mechanics it is 
natural to adopt a material or Lagrangian description of motion in which the 
motion of material particles is traced relative to a fixed reference con- 
figuration. Thus, one can imagine an actual mass of material, the particles 
of which are identified (labeled) in some way. and then one proceeds to de- 
scribe the motion of this mass by giving the spatial positions of each par- 
ticle relative to a specified (generally fixed) frame of reference at each 
time. t. Some of the earlier analyses of special classes of fluid-structure 
interaction problems employed such Lagrangian descriptions, and we mention 
in this regard the 1980 publication by Kennedy and Belytschko (Ref. 14). 

On the other hand, theoretical fluid mechanics traditionally employs a 
spatial or Eulerian description of motion in which the motion of the fluid 
through fixed positions in space is characterized as a function of time. 

Then different fluid (material) particles may occupy the same place in space 
at different times, and the object is to develop the klnematical description 
of motion in terms of these places rather than In terms of the particles. 
Perhaps most of the computational procedures in use for hydrodynamics prob- 
lems employ an Eulerian description of motion, and some of these have been 
applied to problems of fluld-scructure interaction. See, for example, Chang 
and Wang (Ref. 15), Harlow and Amsden (Ref. 16), Wang (Ref. 17), Belytschko 
(Ref. 1), and Dianes, Hirt, and Stein (Ref. 18). 

It is clear that in a general fluid-structure interaction problem, 
neither the Lagrangian/matcrlal nor the Eulerlan/spatial descriptions are 
completely satisfactory. It would be fruitless, for example, to attempt to 
trace the motion of fluid particles in most complex flow phenoo^na (e.g., 
stirring of fluid in containers); moreover, the velocity of the fluid at 
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fixed points in its doaeln is generally the quantity of intereat, not the 
displacement of a particle relative to a fixed point. On the other hand» 
the motion of a solid through a fluid is most naturally characterized using 
a material description, but it is this very motion that alters the spatial 
domain of the fluid with time. 

There are also computational advantages and disadvantages Inherent in 
each of these classical descriptions of motion. In the material descrip- 
tion, the finite element or finite-difference mesh is Imprinted on the ma- 
terial. Thus, with large deformations of the structure, severe distortions 
of the mesh frequently occur, and this has an adverse effect on the numer- 
ical stability, efficiency, and accuracy of most computational procedures. 
This mesh distortion can be somewhat compensated for by using rezonlng tech- 
niques wherein new meshes are drawn on *:he deformed configurations at var- 
ious time Intervals; but these procedures are expensive, difficult to 
Implement, and not completely effective in many situations. The use of an 
Eulerlan scheme to trace both the motion of the fluid and the <^ulld is also 
Imperfect: one must locate material particles of the structure in an 

Eulerlan mesh, and at any particular time the material surfaces of the solid 
will not, in general, coincide with the spatial grid lines defining the 
mesh. Some analysts have, nevertheless, attempted to model the geometric 
changes in Ihe structure with time in an Eulerlan description by using a 
very complex catalogue of material orientations possible in each grid cell 
(see e.g., Chang and Wang (Ref. 13)). The complexity of such procedures, 
and of their implementation has discouraged their use in most fluid- 
structure analysis procedures. One might also mention the presence of con- 
vective terms in the momentum equations for Eulerlan descriptions of motion. 
These destroy symmetry in the resulting stiffness equations and lead to many 
notorious numerical difficulties. While such terms are unavoidable in an 
Eulerlan description of nonuniform fluid flow, their u in the description 
of the motion of solid bodies can lead to 111-condltlonlng of the systems of 
equations governing the discrete model. 
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2.3 MIXED EULERIAN LAGRANGIAN DESCRIPTIONS 

In view of the difficulties noted above, several Investigators have at** 
tempted to develop nixed lagrangian Eulerian descriptions of motion , which 
will be studied in soue detail in the Appendix A. These descriptions gen- 
erally employ, in addition to the spatial and the material frames of refer- 
ence, a referential or grid system that illows one to displace the finite 
element mesh so that it is either fixed, in space, moves with the body, or 
assumes a motion intermediate to these extremes. 

The use of a so-called referential frame, distinct from the material 
and spatial frames of reference, to describe the motion of a continuum can 
be found in several sources on continuum mechanics. A brief discusulon of 
such systems is given by Truesdell (Ref. 19) in his "Mechanical Foundations 
of Elasticity and Fluid Mechanics." However, the Intent of such develop- 
ments does not seem to be to provide a basis for studying the Interaction of 
fluids and solids, nor can one find discussions of kinematics of continua 
sufficiently general to apply directly to interaction problems in any of the 
standard references on continuum mechanics. Interest in "mixed Eulerian- 
iagranglan" descriptions seems to have originated in the computational 
mechanics literature. 

The first attempt at developing computational procedures which employed 
a mixed Lagranglan-Eulerian description appears to have been in the 1964 
papers of Frank and Lazurus ^Ref. 20), and Noh (Fef. 21). These authors 
developed a finite difference scheme for compressible fluid flow in which 
the motion of the fluid relative to an arbitrary moving grid appears in the 
governing equatlcns of motion. These formulations attempt to provide for 
the proper handling of boundary nodes on the fluid-structure Interface while 
allowing nodes interior to the Eulerian mesh to remain fixed and undlstorted 
by the motion of the fluid. Since the resulting formulation ratains many 
features of the &ilerian schemes (e.g. convection like terms), the term 
quasi-Euler^an is also used to describe them. Another quasi-Eulerian finite 
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dlf ferine* Mthod wa« propot«d In 196S by Trullo (kef. 22) with regard to 
the AFTOH hydrodynamic codei and, a decade after the paper of Frank and 
Lasarua (Ref. 20), Aaeden and Hirt (Ref. 23), of Loa Alamoa Laboratories, 
introduced the ALE-technique : Arbitrary Lagrangian-Eulerlan (ALE) scheme, 

which was a finite-difference procedure designed to handle Eulerian and 
Lagrangian descriptions of motion simultaneous!^. More recently, finite 
element codes based on certain aspects of the ALE-strategy have been dis~ 
cussed by Belytschko and Kennedy (Refs. 6 and 7), Belytschko, Kennedy, and 
Schoeberle (Ref. 24), Donea et al (Kefs. 8 and 0 ), Kennedy and Belytschko 
(Refs. 3 through S and 24), Hughes et al (Ref. 23), Liu (Ref. 13), and Liu 
and Ma (Ref. 26). These mlxed/quasl-Eulerian schemes are not without short 
comings: while they provide for flexibility in descriptions of kinematics 

and physics, they involve certain features which lead to the necessity of 
nonconforming finite element methods (see Hughes, Liu, and Zimmerman (Ref. 
25)) and the effects of these bullt*-ln discontinuities on the accuracy and 
stability of finite element calculations is, as yet, not known. 

In the next section, we shall examine the question of appropriate de- 
scriptions of motion in more detail. 
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3. KINEMATICS OF MOTION AND DEFORMATION 


3.1 PRELIMIMARIES 

In aodnrn conclnuua aechanics, the study oi! klaeutlce of continue gen- 
erally begins with a oatheostical characterisation of a continuous body: a 

body "B" is a differentiable aanifold the elements of which arr called par- 
ticles; there is assigned to B a o-finite Borel asasure, called aass of tiw 
body. Thus, B is a natural model of a given quantity of natter as a **con- 
tinuua.” 

Kineaatics aims at describing the notion of B as a function of tine. 

For this purpose, we introduce a tiae scale S c X: and measure the notion of 

B beginning with a fixed instant t > 0 and over a tine xutarval i e S - 

[0,t]. At each t, the particles of B are in one-to-one correspondence with 

3 

points in regions of three-diaensional Euclidian space E : X ; indeed, 
the notion of a body implies the existence of bijective (indeed, diffeo- 

— 3 — 

aorphic) naps k : B Q C H for each tiae where C is the closure of an 

“t 3 — 

open region C X . The regions (or, technically, are called the 

configurations of the body. 

To give aeanlng to the naps and to effect a labeling of the par- 
ticles of B, a particulsr configuration called the reference configu- 

ration, is selected. Typically, fie " i*«*t the reference configura- 

— 2 a o 

tlon is the actual region in X occupied by the body at x 0, but this 
is not a necessary choice of the reference configuration. We proceed to 

3 

introduce in a fixed coordinate on X with origin 0 e and we 
denote by X the position vector of points in 0^. In particular, 1^ 

- and K - Kq, we set 

X - <(X); <: B H. X^ (3.1) 
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and eay that the particle X occupiea the position X la the reference config- 
urer; too of the body. If X^, k • 1,2,3, denote the conponenta of X rela- 
tive to a<»e fixed baaia, then (X^) - (X^, X^, X^) are referred to 
aa the aaterial coordlnatea of the particle X. 

Phyaically, the aituation can be vie«fcd aa thia: He vlah to nane 

(label) the elenenta of B, nei^borhooda of which can be regarded aa actual, 
phyaical plecea of natter. To do thla, we pick one of the reglona In apace 
occupied by the body during ita notion, called the reference configuration, 
and eatabliah in thla region a fixed coordinate ayaten X ~ (^). Nat- 
urally, aince we uaually trace the notion of the body fr<M an initial tine 
X ■ 0, thia reference configuration la ordinarily the actual region occupied 
by the body at x « 0, If X la the poaltlon (in apace) occupied by the par- 
ticle X at X -0, then the correapondence X ■ k(X) effectively aaaigna the 

- 1 2 3 

nunbera (X , X , X ) aa labels (naterlal coordinates) to the particle 
X. Thus, X has the sane label (X^, X^, X^) for all tines x ^ 0. Since 
K la an isonetrlc isororphlsm (relative to the usual Euclide%n netrlc), it 
is usually unnecessary co distinguish between X and its label X in all sub- 
sequent descriptions of the motion of the body B. When the equations of 
notion of B are written in terns of the naterlal coordinates X, we obtain a 
naterlal description or Lagrangian description of motion, as will be 
further expanded upon below. 


Some argue that this is a misnomer since the equations of motion in terms of 
the material coordinates were first given by D'Alembert and not Lagrange. 

See Truesdell (Ref. 19). However, reference to Lagrange here may be due to 
the analogy of this strategy with that employed by Lagrange in his Mechan- 
ique Analytic where he labeled collections of discrete particles and traced 
their motion relative to a fixed spatial frame of reference. This is essen- 
tially w> 't is done here with one fundamental exception: a discrete system 
has a countable number of particles; thus, natural numbers n N can be used 
as particle labels. The Body, B, being a continuum , is nondenumerable; 
thus, a labeling > leme such as the use of triples (X^, X^, X^) of real 
ntabera is needeo to label the material particles. 
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3 

Now, In Addition to the aeterial coordinate frave, points in H are 
Identified by the position vectors x denotes the spatial position occupied 
by particle Ze B In a given configuration 0^, then there wist exist a 
bijective aap k^: B such that 

X - K^(X) (3.2) 

The coordinates (x ) of x relative f.o a given basis are called the 
spatial coordinates of the particle X at tlttc, t. For each tiae reS, there 
exists a unique configuration {2^ of the body, and the family of con- 
figurations, dependent on the real paraaeter, t. Is called the notion of the 
body. Instead of Eq. (3.1) ve describe the motion by a aap <: 

3 

B X [0,t] K of the fora. 


X - v(X,t) (3.3) 

Then the curve x(t) " v(X,t) for finding X, is the path followed by the par- 
ticle X during the motion of the body. Recalling that X ” k(X) and that k 
I s bijective, ve can also describe the notion In terms of the material 
coordinates X: 


or 

where 


X - v(X,t) - v(k"^ (p.t) 

X - x(X,t) 

'' -1 
X ■ O K 


(3.4) 

(3.5) 


Thus, the equation x X(X,t) describes the motion of the body relative to the 
reference configuration In terms of the vectors X or, equivalently, th^ ma- 
terlal coordinates X . When the equations of motion of the body B are 
written in terms of the spatial positions, x, we obtain a spatial or Eulerlan 
description of the motion. When no confusion Is likely, we shall refer to X 
and X Interchangeably as "a material particle." 
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The velocity of the Mterlel particle X at tine t (relative to 0) la 
the vector 


V 


3x (X,t) 
3t 


(3.6) 


If f la a acalar-valuea function of partlclea XeB» the rate at which f 
changea In tlae for fixed X la called the aaterlal rate of change of f . We 
can deacribe thla rate ualng the material tlae derivative 


Dt 


3t 


X fixed 


(3.7) 


Thua, 


Dx/Dt 


(3.8) 


If g la given aa a function of the apatlal coordlnatea x X(X,t) of 
partlclea X, then 


^ ^ 38(x.t) 3X,^ 

Dt at 3Xj^ 3t 


(3.9) 


where X|^ denote, for Inatance, coccdinatea of x relative to a fixed baala 

In R (aay Cartealan) and repeated Indlcea are auamed, k « 1,2 N. 

Thua, If Vj^ are the correapondlng components of y. 




98 4. V 98 
3t ^ ""k axj^ 


(3.10) 


3.2 THE QUASI-EULERIAN DESCRIPTION 


We now set out to derive a description of motion which la sufficiently 
general to encompass both the material and spatial descriptions as well as 
Intermediate mixed descriptions that may be appropriate for fluid-structure 
Interaction problems. We continue to employ the notations and conventions 
Introduced earlier: for a material body, B, a family of smooth bljectlve 

maps such that 
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II 

Kq! B -V C R * Is the reference configuration 

X ■ !5q(^) “ labels of material particles 

■ positions of particles In the reference configurations 

N — 

B -> ( 2 ^ R , Is configuration of the body at time t. 

X - !St (X) - (|c"^ (X)) X(X,t) 

- the spatial position of particle X at tine t 
N - 1, . or 3, 

In addition to these quantities, we Introduce a smooth bljectlve map ^ from 
N 

X (0,T) into S , such that Is the Image of 41 at time t: 

For teI0,T], H and we write 


X - f(y»t) (3.11) 

or 

y - 4 )“^ (x,t), X e t e [0,T] (3.12) 

Since the map X Is invertible, we can also write 

y - (X (X,t), t) = 4j(X,t) (3.13) 

Ihe vectors y are said to specify referential positions of particles at 
various times, t. y refers to the position x of particle X at time t relative 
to a moving frame of reference in $2 with origin 0. Thus, we Imagine 0 


H 

In this respect, we depart from certain mixed Lagranglan-Eulerlan descrip- 
tions found In the literature which hold i 2 fixed for all times but allow 
n(*S2_) and n to be time dependent; see, for example, Hughes et I 
(Ref. 25), Liu (Ref. 13), and Liu and Ma (Ref. 26). 
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moving with time through space and y ■■ ^(x»t) a vector from 0 to the spatial 

- - ~ k k k 

positions X of particles. We also Introduce coordinates X , x , y » k 
"I. N, relative to fixed set of basic vectors In each of the respec- 

tive d(»alns. These notations and the geometrical situation are depicted in 
Fig. 3-1 for the case In which o and 0 coincide and are fixed In space. 

The coordinates y are frequently called the grid- or mesh- coordlnates , 
for reasons which will become clearer later. It Is worth noting here, how- 

A 

ever, that if we lay a f** 'ed mesh on n, the choice y=X (l.e., rj/ I the 
Identity) yields a Lagranglan grid whereas the choice y ■ x produces an 
Eulerlan grid. 


In addition to the particle velocity field, (Eq. (3. A)), we Introduce 

g 

the grid velocity v , which Is the rate at which a "grid point" y moves 
from a fixed position x In space; l.e.. 


X - (3. 14) 

y"fixed 

The so-called difference velocity v^, defined by, 



v 


D 



(3.15) 


is then the velocity of a material particle relative to the moving grid n. 


As preparation for a major transformation rule pertaining to material 
derivatives, we establish the following lemma: 


Lemma 3.1, Subject to the conventions and assumptions established 
above . 
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X = X(X.t) 



Fig. 3-1 - Various Regions and Coordinate Frames Characterizing 
the Motion of a Material Body, B 
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3t " ' 9Xj • 


(3.16) 


where i|» ^ Is the Inverse of the nap ^ of Eq. (3.11) and the sunnatlon 
convention Is used (j, k - 1, 2, ... , N). 


Proof ; 
Eq. (3.13), 


Let denote Cartesian coordinates of y. Then, according to 


Differentiating this expression with respect to tine holding y|^ fixed 
yields 


as asserted. 


JT ” 


0 


^^k^ j 9’^'k^ 

3Xj * t 3t 


I^t y be scalar field given as a differentiable real-values function f 
of the material coordinates X and time t: 

V - f(X,t) (3.17) 

Then the material rate of change of p is the rate p changes in time for 
fixed X: 


Dp . 

Dt " 3t 

The following result allows us to compute the material time derivative of 
In terms of the grid coordinates and the difference velocity. 

Theorem 3.1. Let the conventions established earlier hold and let p be 
given by Eq. (3.17). Let 
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g(y,t) - (y»t), t) 


“1 


h(x,t) - fQl; (x,t), t) 


0.18) 


lh«n 


5E 

Dt 


3g(y,t) 3h(x,t) « 

-3T- + -Tij- 


(3.19) 


Proof; The Important Idea to be kept in mind here la that we wish to 
write y as a function of y but compute its time-'rate-of'- change holding X 
fixed . Wa have , 


Du 

Dt 


8(y,t) 3g(^”^ (X(X 


9t 


ayi 


,t).t).t) 

y 3xj at at ) 


But 


9h(x,t) 9h((j>(y,t),t) 9<^”^ 9g(y,t) 3<|)“^ 


9x 


j 


9y^ * axj 


3y^ ’ 3x. 


and, from Lemma 1, 


3g(y,t) 9(|.j^^ 

3yu * 


-1 


3g(y,t) 9(()j^ g 9h(x,t) 


9yjj * 3Xj j 


9x, 




Thus, 


Du 

Dt 


3g(y,t) 3h(x,t) 




( Vj-vp 


as asserted. 

It is Important to note that Eq. (3.16) reduces to conventional Lagran- 

glan or Eulerlan descriptions with appropriate choices of the coordinates 

k . . , 

y or the map 


16 


LOCKHEED-HUNTSVILLE RESEARCH & ENGINEERING CENTER 


UlSC-HBSC TR D66728S 



ORIGINAL PA2E IS 
OF POOR QUALITY 

Lagranglant 

♦ - X, V - — - y. V - 0 


Dy 3g(«(X.t),t) 3f(X.t) 

TJF at at 

Eulerlan: 

♦ - I, V® - Xl x - 0. V® - V 

‘ ^ ®‘|x fixed 


jjy 3g(;»t) 3h(x,t) 

Dt at ^ ax. ^1 


3,3 OTHER KINEHAXICAL EQUATIONS 


Some isiapllficatlons In Che developments can be realized by considering 
the referential coordinates to coincide with the material coordinates of a 
particle X at time t * 0 and to regard the "grid" as moving relative to the 

Q 

reference configuration 0 at an arbitrary velocity v which Is not di- 
rectly dependent on the particles. If X^, y|^ are Cartesian compo- 

nents of X, X, and y relative to fixed bases, we have 




G ^^k 
^ (X.t), \ 


Let 


J(X,t) 


det 


3X“ 


Then we have the following results: 


Lemma 3.2 


H 

Dt 



(3.20) 


(3.21) 


(3.22) 
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Proof: Lot 




and lat danota th« altamatlng tanaort. Than 

J - datfy, 1 - r 7 7 7 

F ® ® ^1,0 ^j,a ^k,Y 


and 


3J 


1 _aee _ 
Ty;;;) - 7= ^i.1r 


since 


’^k 3l“ 


3x“ ■ 


we have 


3JL . “*««<»■• lyt.,1 _ 1 _ p6a 


3y 


j 


d«t Iyk,a^ ^ 


Thus, froa Eqs. (3.23) and (3.24) 

3J 


^<='r,6> 


dx 


e 


- cofactor y^^^ “ 


Also note that 


A - ^ !!l . G 

at ^r,9 " * at " ^r,e 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


Thus, 
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3J _ 3J ^^r.e 

^ ■RSvTT- “§t- 




3,3. Given a differentiable real~valuea function g ■ g(y,t) ■ 
g(^(X»t),t), we have 


( 1 ) 44 ^ ^ 8 ^ G 

^ ^ « 3t 3y. % 


(3.27) 


Slwllarly, 


(2) Sljt-li . J(|4 + Y . (gyC)J 


Proof; Condition (1) is obvious. To obtain (2), note that 

-Gv _ a 


I . (gv ) £ ^ (g(y.t) v^(y,t)) 


(3.28) 


G G 

g 7 . V + V . 7 g 


JV . (gv“) 


C G 

gjy .V + Jv . 7g 


g |7 + Jv®.7g (3.29) 


where we have aade use of Leaasa 3.2. H?us, froa Eqs. (3.27) and 3.29) 






- Jlf + J 

f 

- J -»■ J7 . (gv®) 


Is asserted. 
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W« vlll ut« that* results la ths osxt ssctlon to dsrlvs s^ustlons of 
■otlon In ths grid (rsfsrsntlsl) coordlnstss y^. 

3.4 A SPECIAL REFBISIITIAL DEFCKIFTION OF MOTION 


A fslrly gsnsrsl discussion of klnsnstlc^ in s grid or rsfsrsntlsl sys 
tsa wss given ssrllsr. Ws will now focus on spsclsl extensions end applies 
tlons of those Idsss which lend to s convenient freaework for treating 
practical prohleas in fluid-structure Interaction. The structure of this 
kineaatlcal description is outlined In the following five steps: 


1. Spatial Ftmw. Ue establish an absolutely fixed (spatial) 
Reference fraae in ths position vectors x are defined 
by their Cartesian coaponsnts x|(, ^ h ^ N. 

2. Material Fraae. At tlae t ■ 0, a aaterlal body B occupies a 
region S’ C 1^** we use as labels of the particles of B 
the coordinates, of their positions In this refsrenco 
configuration; the corresponding places of particles 

are Identified by vectors X. 

Motion. The notion of B is, as usual, defined by the specifi- 
cation of ths position X e of each particle X at each 
tlae, t,0^tjf^T: 

X - X(X,t) (3.30) 

The action ^ Is assuaed to be a differentiable bljectlve aap 
froB n into n Into C R*^ at each tlae, t. 

4. Crid Positions. We introduce an arbitrary, differentiable, 

Injective function 0 x |0,T] R^ x [0,T], such that for 
each { e [0,T], the range of ^ Is a region 0^ C The 

"position" vectors 

y - <KJ,t), y e fit (3.31) 

are said to define the grid positions In R^. Note that y 
Is a position vector of a point (place) in R^. These posi- 
tions depend only indirectly on the locations x of aaterlal 
partleals, vlt, 

y - 4(X”^ (x.t),t) - ^(x,t) (3.32) 

5. Dlsplaceasnts and Velocitieo. One can define the particle 
displaceaent u and the particle velocity, y> hy 
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u(X,t) - X - t • X(X,t) - X 
Dx DX<X,t) 


(3.33) 


- dT - 


dc 


•nd th« grid dltplaccMnt and grid vtlocity by 
uC(x.t) - y - X - ♦(X.t) - X 

^ ^ ^ m m ^ -- 




3(^(X,t) ay 

* ^ ■ «J1 

3t n 


(3.3A) 


X fixed J 


Iheae quaatltlea and notatlona are llluatratad In Fig. 3-1. 

3.3 JACOBIANS AND TIME RATES 

Let denote coaponente of y and X relative to a fixed baala. We 

denote by j the Jacobian of the tranaforaatlon x|-^ A(R«t) for each t: 


J(X.t) - det 




ax 


det 


^Vx.t) 


9X 


(3.35a) 


Llkcvlae, we denote 


J(X,t) - det 


9X 


■ det 


94^(X.t) 


9X 


(3.35b) 


Then, as was proved In Section 3.3, 



(3.36) 
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Note that th« operator V la the saao spatial gradient operator la both of 
thes expressions owing to our definition of In both expressions, the 
tlBe-rate**of -change Is taken with X fixed. 

Next, w let g denote values of a given real-valued function of ;<ar- 
tlcles and tlae, and Introduce the notations 


where 


g ■ value of field at particle X at tine t 
- g(X,t) - g(x,t) - g(y,t) 


g(?»t) 

g(y.t) 


g(x"^ (x,t),t) 
(y,t),t) 






( 3 . 3 ?) 


< 3 . 3 ^) 




Dt 


aaterlal tlne-derlvatlve 
tlae-rate-of-change of g 

ii 




If + y . vi 


• vg 


3t 


of g 

for fixed particle X 


( 3 . 39 ) 


where it is undereto<)d that 

9g _ the tlne-rate-of-changc of g holding the 
3t spatial position y fixed. 


( 3 . 40 ) 
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Thua, while the fons of the functiona g and g aay be different! the par- 
tiala 3^3t and 3g/3t both represent the re«.e g changes in t at a fixed 
spatial position. 

Coablnlng the above results, we arrive at the equations 



(3.41) 


Also note that 


(3.42) 


We now make a fundamental observation : the partial derivatives 3g/ 

3t and 3g/3t represent tlme-rates*-of-change at fixed points In space. The 
values of these rates coincide If we take ^ « X and v - 0. Thus, 


3(Jg) 

3t 


X fixed 


+ Y • 


,G, ^ 


-^■^1 - J(|f- + V® . yg + gV . V 

lx fixed ' 


This observation was apparently first made by Donea et al (ftef . 10); see 
also Donea (Ref. 12) and Belytschko and Kennedy (Ref. 27). 
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38(X(X,t),t) 

3g($(X.t).t) 

3t 

-el “ 

y-0 


(3,43) 


This relationship proves to be crucial in deriving local equations of notion 
In referential coordinates since it enables us to transfom any standard 
Bulerlan fora (spatial tiae rates) into a corresponding tlae derivative In 
grid coordinates. We will exploit this Idea in Section 4. 
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4.1 INTRODUCTORY REMARKS 

We consider the motion of a material body subjected to body forces of 
intensity, b, per unit mass and surface tractions, S, on a portion, 
of its boundary. We focus our attention on control volumes H C IR with 
boundary 3J2 » U the velocities being prescribed on 3fij^ and 

the tractions on 9^2* Alternatively, to obtain equations of motion In 

A 

referential coordinates, we consider a control volume, moving 

^ ^ ^ ■ I . ^ . ^ 11 1- 

with the grid velocity, v , and with boundaries, ■ 3Jij^ U 3^2* Com- 
ponents of vectors are referred to a fixed orthonormal basis characterizing 
the spatial coordinates, x^^. Indlclal notation and the sumnatlon conven- 
tion are used in some of the relationships which follow. The following 
additional notations are Introduced: 

p ■ the value of the mass density of the body in the cur- 
rent configuration, 

0 - p(x,t) - p(x(X,t),t) - P(X,t) 

0 a the Cauchy stress tensor, with Cartesian components 
relative to spatial coordinate directions of 

oij “ oij(x,t) - Oij(X(x,t),t) - Oij(X,t) 

b - the body force vector per unit ma&s in the current 
configuration with Cartesian components 

bi a bi(»,t) a bi(X(X,t),t) a bi(X,t) 

e a the Internal energy per unit mass In the current con- 
figuration of the body, defined by functions 

e a £(X,t) a e(X(X,t),t) a F^.t) 

Other notations will be Introduced later. 
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4.2 LOCAL EULERIAN FORMS 

The local Eulerlan (spatial) forms of the equations of motion of an 
arbitrary continuum are: 

1. Conservation of Mass 


+ 7 . (py) - 0 


(4.1) 


2. Balance of Linear Momentum 


apv. 






8o 


ki 


(4.2) 


3, Balance of Angular Momentum 


'ij - “"jl 


4. Conservation of Energy 






(4.3) 


(4.4) 


where 


total energy per unit mass 
1 


j : 


V . V + E 


(4.5) 


heat working 


9x. 


+ pr 


(4.6) 


and are the components of the heat flux vector and r is the heat sup- 
plied per unit mass per unit time in the current configuration. For sim- 
plicity, we will take Q ■ Q In subsequent developments; but it should be 
clear that the addition of Q and thermal effects produced no significant 
complications in any of the following developments. 
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4.3 local referential FORMS 


A direct application of relations (3.41) and (3.43) derived in the pre 
vioua section and the above Bulerlan forms leads to local eqxiatlons of mo- 
tion In the grid-coordinates, Y|^. 


The superposed caret ('') denotes functions of the referential coordi- 
nates and time. Thus, for example, the mass density is given by 


p - p(yft) 

where It Is understood that 

P(y,t) - p(KX,t),t) 

= p(X,t) = p(x,t) 


etc. , 


(4.7) 


1. Conservation of Mass. In view of Eq. (3.41a), 


ajp 

at 


X 



where Indicates that the time derivative Is taken with X held fixed. 
But, according to Eqs. (3.43) and (4.1), 


Thus, 


9p(y.t) 

at 


3p (X» t ) , t ) 


at 


V ■0 



A 

Ik 

at 


X 


jv . PV® 


(4.8) 
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where v is, again the difference velocity 


D G 

V*' c V - V 


(A. 9) 


Alternatively, since 




- PjV . V® + j 1^ 


ve have 


^ I 

i£. 
at , 


- pV . V® + V® 1^ + pV . (y® - y) 


or 


ap. 

at 


A A 


V . 7p - pV . V 


(4.10) 


2. Balance of Linear Momentum. Let 


pb. + 


3 a 


kl 


1 ayi 


(4.11) 


Then a calculation similar to that leading to Eq. (4.8) yields the momentum 
equation in referential form. 


JpjVj^ 


at 


:d. 


J ^p^k ^i> JPk 


(4.12) 


Similarly, by expanding the left-hand side of this equation and using 
Eq. (3.41), we have the equivalent equations. 
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n 


p^k^^Pi “^1 


- :) 4: 


(4.13) 


According to Bq. (4.8), the term In brackets vanishes if mass Is con- 
stant locally. 


3. Balance of Angular Momentum. Denoting 


Oij(y,t) 


oij((j>(X,t),t) 

Oij(X,t) 

Oij(X“^(x.t),t) 

cfij(x,t) 


angular momentum Is balanced locally If 


(4.14) 


ojj(y,t) - Sji(y.t) 


(^.15) 


4. Conservation of Energy. By following the Identical process used to 
obtain Eqs. (4.8) and (4.12), we arrive at the following referential form of 
the conservation of energy. 


/N 

9jpE 

9t 


X 




OEV®) 


- JloVk •* -9y- 


+ 3^ (PEV®)) 


(4.16) 


Since E ■ 1/2 v^ + e (v^ ■ v . v), we can ;ilso write this result In the 
form 
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?jjki 


dt 


9V. 


- J 


•5~ (P.vp + ^ 


V 

2 


“ jv. 


4 


3pvJ 
X ' ^ X' 


3t 


X ^ J 3yj ay 


Ik 


(4.17) 


If mass Is conserved, the term with a single wavy underline vanishes, 
by virtue of Eq. (4.8). Likewise, the term with double wavy underlines re- 
duces then to the local momentum equation (Eq. 4.13) which also vanishes. 
Then If linear momentum is balanced, we have 


a(jpe) 


3t 


av,. 


X 


J 4 (Pev”) 


(4.18) 


Finally, since 


/N 

ajpe 


at 


- J 


, ^ 


3t 


+ pejv • V 


we have 


5 >Ji + V® - de — 
3t „ °kl 3y. k 3y, ^ 3y, 


(4.19) 


5. Summary. In summary, the local equations of motion In referential 
form are given by 
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(4.20) 


wherein it is understood that the tlne-derlvatlves on the left side of the 
equality are computed holding X fixed and that the quantities appearing on 
the right side are regarded as functions of the grid coordinates, 


4.4 GLOBAL FOBMS 


Let dv be a differential volume element in 0 and ds an element of 

surface area of the boundary, 90. , with unit exterior normal, n. The con- 

^ ^ ** G 

trol referential volu^iie, 0^, Is moving with the grid velocity, v , rela- 
tive to the fixed spatial frame of reference as before. Let dv^ denote a 
material volume element In the reference configuration, so that 


dv - I det ^ ]du^ - jdv^ 
a 


Finally, let G be a quantity to be conserved in a physical process and sup 
posed G Is given by 


G(t) 


/ 


g(y,t) dv 


Then the time-rate-of -change of G is 
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at 


3 (y,t) dv 




/ 


Ml dv 

at o 


f (|f ) + 7 . v®g)J 

y* dv + y ^ 


Hence, an application of the divergence theorem yields 


dG(t) 

-~ST^ 


/ If ./ Y® . ! * « 

'' an*' 


(4.21) 


1. Mass. Since 


at 


- V . pv 


J ‘ - J 


3i' 


Hence, if M(t) is the total mass of at time, t. 


dM(t) ^ d 
1[t~ 


A A 


Jt ■ 


/ 


TV A A 

py . n ds 


an. 


(4.22) 


(4.23) 


2. Linear Momentum. If P(t) is the total linear momentum of at 
time, t, then an application of Eq. (4.21) yields for the global form of the 
rate- 0 f-charige of P, 
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dP(t) 

IT- 


d 

art 


/ A A A 

pv dv 


an 


J * • §) 




/ 


(dlvo + pb) dv 


(«.24) 


J^!££FX^ Lik«wi*e, the rate-of-change of total energy le 


. d 

dt IF 


, /* pE dv - r pb . V dv + /* (p 

^ anr 


(pE v" . n 


A A A 

+8 . y) ds 


(4.25) 
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5. VARIATIONAL PRINCIPLES 


5.1 SPACES OF ADMISSIBLE FIELDS 


We ehall introduce ebetrect linear epacea of adalasible densltlee» 
velocitlea, and Internal energy: Let H denote a "control voltue" in the 

n 

referential ayatea aovlng with a velocity field, v , relative to the fixed 
apatlal fraae and let 30 denote Ita boundary. The boundary, 30, ia further 

A ^ C 

decMpoaed into portiona, 302 where the velocltlea y, v and 

A A A ^ 

the tractlona aj^ ■ ^ki^l preacrlbed, reapectively, n being a unit 
noraal to 30. We have thua, for each t e [0»T], 0 ■■ (4>»t) 

lot 


V ■ apace of adaiaaible velocltlea 


r - - 

J ^ dv 

0 


< 00 , V ■ 0 on 30 


R ■ apace of adaiaaible denaitlea 


(v^ . ^ I ^ 1®»T] 


E - apace of adaiaaible Internal energy denaitlea 


/ 3v, _ 

('/'Ori . y + ij/y , dv| < », 


V^CR, Vv e V, Vt E 10,T] 


(5.1) 


(5.2) 


(5.3) 
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aethemetlcel eignlflcence. ttore specific properties of these spaces cen be 
established once constitutive equations for the aaterlal have been defined, 
typically, the spaces V, R, and E will be Orllcs^Sobolev spaces such as 
R ■ W**^ (n, L^(0,T), etc. We shall give acre concrete definitions of 
these spaces In later developaents. 

5.2 VARIATIONAL PROBLEM IN REFERENTIAL FORM 

We now couslder the followlug variational problem: 



It la easily shown that any solution of the equations of motion (4.20) 
Is also a solution of the variational Eqs. (5.4). Conversely, any suffi- 
ciently smooth solution of Eq. (5.4) is at least a weak solution of 
Eq. (4.20). Thus, Eqs. (5.4) represent a set of variational equations 
equivalent to the equatlonri of motion Ic referential form. 


35 


LOCKHEED-HUNTSVILLE RESEARCH 6 ENGINEERING CENTER 



LKSC-HREC TR D867285 


ORIGINAL FACZ t:? 

OF POOR quality 


6. FINITE BLEMENT MODELS 

6.1 THE DI8CBETE VARIATIONAL PROBLEM 

The construction of finlto elcnont aodels of th« genornl fluid- 
•tructura intaractlon probloa aabodlad in Eqa, (5.4) follows tha standard 

staps: 

A 

1. Tha doaain^ 0, is partitionad into E finita alaaants such that for 
aach t c [0,T] 

n ■ ^ a ■ f 

a • * * 

2. Piacawisa pol 3 moBial shspr functions sra dafined over each alawant 
which rpovide a basis for local approxiaationi of p , u, and e. Typically, 
these have tha properties 

a 

"h 


where R » N , and K denote the nuabers of element -degrees-of -freedom 

* ® ® N M L 

for tha respective local approxiaations; p^, and denote nodal 


R_ 


T, 

h«l 


p*(t) sj(y) 


E 

M-l 


u^(t) 4>*(y) 


( 6 . 1 ) 


E 

L-1 


E^(t) ipj(y) 
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6 0 6 

values of p , u , and f at tlae t, and the local shape functions 
h 'h h 

usually satisfy conditions of the type 




.6/ N. 




6^, 1 ^ M, N 

5 “. 1 M, N < Ng \ 

5 “. 1 < L, N 


( 6 . 2 ) 


N H N 

with y , y^, and y^ the element coordinates of nodes corresponding to the 
local approximations of p, u, and c, respectively. 

3. The shape functions are designed so that they match at inter- 
element boundaries so as to produce global basis functions 


Bj(y). $j(y), 4»j^(y) (1 i < R. 1 1 j 1 N. 1 £ k ^ K), 

which are defined over the entire domain n and which provide basic functions 
for flnlte-dlmenslonal spaces R^, and respectively, and so that 


B. 






'N 


(6.3) 


In conforming finite element approximations, frequently have 


C R, C V, C E 


(6.4) 


where R, V, and E are the spaces of admissible functions introduced in the 
previous secton. We shall not, however, restrict our analyses to conforming 
finite elements. 
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6.2 SEMI-DISCRETE MODEL 

The seal -discrete, Galerkln, finite eleaent approxiaation of the varis' 
tional problea (Eq. (5.4)) is characterized by the following discrete prob- 
lea: 


Given f, V , and S, find (p. , u. , e, )e 

_h „h~ _h . , 

R , y , E such that 


/ ^Ph - /* D 1 ‘ 

w * W • ?h> • J *h 


r , =Sh h B 

. J <»h^h • ^ % ''l V 3P— > 


hk 


/ 


3v 


(o, 


hi 


+ Pv K ) dv + 


ki 3y. ^h “1 'hi 


an 


/ 






(6.5) 


r , . ^"h\ ^ , -^"hk D 

J at ’^h Ph ^h 3y. '^h \k dy. 


^\k 

k 


^Ph^h 


au 


hk 


-o, , ) dv - 0 

ki ay^^ ^h' 


where 


V(*J'h. vh, 'i'h) E R^ X X E*' 


D 

^h 


V® - V 

^ ~h* 
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7. INTERFACE CONDITIONS 

The main difficulty in enploying a general ALE-type description of 
notion for fluid-structure interaction problems is the specification of 
boundary conditions at the Interface of Eulerlan and Lagranglan meshes and 
at free boundaries. If V denotes any material surface on which such condi- 
tions are to be imposed, then a necessary condition to be met is 

v^"vonr (7.1) 

n 

G G 

where v ■ v . n, v ■ v * n, and 

^ n n# 

Q 

V - the grid velocity 
*** 

V “ the particle velocity 
n " a unit normal to T 

*v 

To generalize this condition, consider the situation shown in Fig. 7-1 

G 

in which a nodal point N is assigned a given trajectory and velocity v . 

The material surface, F, is assumed to be given by the equation 


If 


1^ ~ unit orthonormal basis vectors 

r “ position vector of points on surface T 


then 


n 


fl + <t>’' 
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r~ Free Material Surface » 



Fig, 7-1 - Geometry of Free Surface 
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oRiar^-:’. ?. 
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■ r. 

t X,’ '-»» 

JALm"< 


and 


V - V 
n n 


h * ''2l2> • ; 


(1 + 


2 TTTI 
(1 + 4 >' ) 


If denotes a coordinate tangent to the node trajectory path, then 
0, and ve have the condition 



2 dx. 


(7.2) 


In general (for three dlnenslons), If r Is given by an equation of the 

form 


Xj - (Xj.Xj) 


the Interface condition Is 


G d}p 

If r Is a fluid-solid Interface, we have 

G Fluid Solid 

V - V ■ V on 

n n n 


(7.3) 


(7.4) 


at all nodes on T. 
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8. A LINEAR PSI-MODEL FOR NEARLY INCOMPRESSIBLE VISCOUS FLOWS 
8.1 GENERAL 


We would like to construct a formulation of the fluld~8tructure Inter- 
action problem which has the following features: 

1. Is applicable to the problem of a viscous Incompressible or 
slightly compressible fluid Interacting with an elastic solid. 

2. Is characterized by linear equations. 

Unfortunately, linear FSI-models for compressible flow altaost exclu- 
sively deal with the small~perturbatlon acoustic approximations of plane or 
spherical waves Impinging on an elastic body and these problems are of se- 
condary Interest here. We shall therefore, consider a model which arises 
from a penalty treatment of the continuity equation for Incompressible vis- 
cous fluids, thereby allowing for a non-zero "bulk viscosity” of the fluid. 

8.2 GOVERNING EQUATIONS OF THE FLUID 

For the general Stokesian fluid, the governing equations of motion are: 


Continuity: 


Momentum: 


f \ ^k P^^k,k ■ 


^,.1 - P^'-^kj,.! 


( 8 . 1 ) 


( 8 . 2 ) 
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Constitutive: 



(- TT 


+ A d ) *5 , . 
rr' kj 


+ 2u d, 


(8.3) 


vhere 

p * the fluid density 

■" Cartesian components of velocity 

bj^ “ body force density 

a, ■ components of the Cauchy stress tensor 

J 

TT “ the thermodynamic pressure 
**<-1 " T ^ deformation rate 

X, y * deletatlonal and shear viscosities. 

The thermodynamic pressure tt Is given In terms of p and the absolute temper- 
ature 9, by the equation of state of the fluid. Since we are Ignoring ther- 
mal effects here, we have 


Equation of State; 


TT - TT (p) 


(8.4) 


For polytropic gas, for example, 

IT (p) - 

where p * constant and y Is a material constant. 

Upon substituting Eq. (8.3) Into Eq, (8.2), we obtain the Navler-Stokes 
equations for Isothermal, viscous, compressible flow: 


P 




Pb^- 


TT + (X+ y) v 


.1,jk 


+ y V, 




( 8 . 6 ) 
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which can also be written In the vector fora 

pv ■ b-VTr+(X+2y)V (V.v)-yV xVxv (8.7) 

where v la the aaterlal tlae derivative of the velocity; 

. Dv 3v 

V - - -a~+v.Vv (8.8) 

8.3 A PENALTY FORMULATION 

Let P denote the vector defined by 

pP « p(v - b) + yV X curl v (8.9) 

Then the equations of motion can he written compactly as 


p ■ ~ dlv V 

pP “ kV dlv V - 7tt (8.10) 

^ ^ ^ 

2 

where k ■ X + y Is a "bulk" viscosity for the fluid. In the case of an 
Incompressible fluid, these equations reduce to 

dlv V ■ 0 

P P - V- TT (8.11) 


where tt Is now the unknown hydrostatic pressure. 


A penalty approximation of the Incompressibility condition (dlv v ■ 0) 
in Eq. (8.11) yields alternatively, the modified momentum equation 


pP 


-1 


V dlv V + kV dlv V 


( 8 . 12 ) 
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vh<>rfi' e Is a positive number and one hopes that 
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with 


dlv V -► 0 
-e 


■t: -► TT 

e 


as £ 


0 


(8.13) 


TTP - 



dlv ve 


(8.14) 


and the convergence In Eq., 


(8.13) Is In (H’(f2))” for v^ and L^(fl) for 


It Is clear that 


1. The penalized momentum (Eq. (8.12)) corresponds to that of a 
compressible fluid with an equation of state given by Eq. 

(8.14). 

2. If this compressible fluid Is characterized as a haratroplc (?) 
gas according to Eq. (8.5), then the continuity equation Is of 
the form 

p « e PoP^+1 (8.15) 

3. Over any djraaln on which p.^ ■ 0 (p varies only with 

time) we have 

o*' p 

that is, the density in can be determined by a quadrature. 

We shall employ the approximation®. (Fq. (8.12) and Eq. (8.16)) In sub- 
sequent discussions. 


8.4 GOVERNING EQUATIONS OF THE SOLID 


The motion of an elastic solid Is governed by the classical equat:lons, 


(E 


Ijkil ^8 ^i 


^ "l 


(8.17) 
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where 


•iM 


“k 

Ps 

fl 


Ok.:. . p., 

OF POOH quality 

the elastic constants of the material (Hooke's tensor) 
exhibit standard sytsnetrles and elllptlclty 

the Cartesian components of the displacement vector 

the mass density of the solid 

the components of body force of the solid 


8.5 BOUNDARY AND INTERFACE CONDITIONS 

Consider the geometrical situations and notations Indicated In Fig. 
8-*l. In this case, the FSI-Problen Is governed by the following equations. 


Fluid 


-1 


pv T (e - k) V div V + yV X curl v ■ pb lnfj_ 


If 


V on r, 


FI 


“ s on r, 


F2 




(8.18) 


Fluid-Solid Interface 


u on Ft 


(8.19) 


Solid 


V.Eu+p f “ Pu onJ^ 




S ~ 8 - 


u “ 0 on r 


si 


n . E u ■ Sg on T ^2 


( 8 . 20 ) 
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f 

where o (v) it defined in Eq. (8.3) and n la a unit vector exterior and 
normal to the surface The Interface condition, Eq. (8.19) charac- 

terizes the ao-called no-slip condition ; l.e., the fluid is to adhere to the 
solid at the interface: to do otherwise would suggeat a discontinuity in 

traction aero ta the interface. 

In certain instances which suggest that this no-sllp condition be 
relaxed, we can use instead of Eq. (8.19) the slip-interface condition 

u - V ; ( j - v) . T ■ . T on F- (8.23) 

n A ^ X 

where u ■ u . n, v ■ v . n, t is a unit vector tangent to F-, u is a 

D ^ X 

film coefficient (which may depend upon temperature and density of the 
fluid) and is given by Eq. (8.22). 

8.6 VARIATIONAL FORMS 

Let Ve (C^(jy_))^ and (/?f (C^(H ))^ be sufficiently smooth 

-w r 8 

vectors defined over?!- and 0 , respectively. Taking the inner product 

r 6 

of the momentum equations of the fluid and the solid with ip and ip , respec- 
tlvely, integrating over the respective domains, and applying the Green- 
Gauss theorem yields: 


f + f N ♦k.i + k. ¥ 

r\ ^ 




I 


o‘,k(v) Iij d. + / Pb . ♦dx 

Sn- ■ 


/ 


(8.24) 


. /" p. \ k 

* n 


IjkJ. k, i,J dx 


f p, f . (?dx + J "i *’k 
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wherein 


k * X + y - e 
c 


-1 


■ "’'l.k + •'c - * 


o..(u) • E, u 

kl kire r,e 


Ik 


(8.25) 


on DQp we teke 


♦k ■ ° ''n ''r ”ki<y> ”k ■ ®" ^T2 


Thus, 


/ “k *i 


dU, 


- f 


n 


F2 


Slmilar.ly, with ■ 0 on ^si ^ 


/ °w<“> \ 1 “• ■ 


3fL 


/ 


®k "^k <>• 


S2 


Thus, we jrrlve at the variatlona.i boundary-initial -value probleD of 
finding v(t), u(t), te[0,TJ, such that 


*k ♦k.i 




/ 


\ dx + / Sj^ da 


/ 


F2 


o7i-(u) 4 dx 


^lk'~' ^k,l 


I 


dx + 


PS ‘k'^k 






S2 


(3.26) 
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for all aufflclently smooth functions ^ which vanish on ^ 

I’gj rsspsctlvsly. 

8.7 FINITE ELEHENT MODELS 

We partition IT ^ ^ Into finite elements In the usual manner. 
The approximate velocities and displacements are of the form 




(8.27) 


where v*^, u'^ denote values of v^, and u^ at a nodal point J in the 
respective fluid and solid meshes. Introducing these approximations Into 
Eq.(8.2S) yields the corresponding element equations of motion In terms of 
the nodal values. The final system of equations for the discrete model is 
of the form 


Jff 2fi 



«IF ”ll ”lS 





0 


2si 5ss 







-F 



S 

h? -II - 


-I 

0 0 0 


• 

“s 

W 1 u 




“ 0 0 0 “ 


1 

b. 
> t 

1 


■* b ” 

0 K,.. K „ 


u.. 


f^ 

^ -IS -SS 




-I 

0 K.., 




f 

- ^I -SS 


-vS 



» — 




_ 


(8.28) 
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Here the matrices 

-^r 


are 


are mass matrices, 1^, 
fluid "stiffness" matrices, a*"* stiffness matrices for the 

solid. The vectors v and u of nodal values of velocity and displacement are 

•V 

partitioned Into column vectors corresponding to nodes on the interior of 
the fluid/solid mesh and nodes on the fluid-structure Interface: 


- Vf- 


-Ui- 


, « * 


v» 


U^ 

_ I _ 


_ S _ 


The no-slip Interface condition enters the formulation by setting Vj ■ 
u^. Hence, the division of nodal-point degrees of freedom In Eq. (8.28) 
corresponds to a convention of the type Indicated In Fig. 8-2. 



Interface 



Fig. 8-? - No-Sllp Interface 
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Since Mpp 
tlons) in Eq. 


Is Invertible, the second equation (the interface equa 
(8.28) can be written 


— «• •• — f — f * * 

M__ u_ + M ii + K V + K + K__ + K__ u_ 

-.II -I -IS ~s IF F -II -I -II -I -SS -S 


!i 


where 


M - M - M M~^ M 

”ll %I -IF -FF ”fI 


^IF “ -IF " -IF ^ ^ 


^II " -II ” -IF Sf ^I 


Cl * Cl “ ?!if ^ - 




(8.30) 


(matrix condensation). 
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9. SPAR MODIFICAriONS FOR INCOMPRESSIBLE VISCOUS FLOW CALCULATIONS 
9.1 ORIENTATION 


In the standard program, SPAR computes the stiffness matrix k for the 
four-node (E41) element by the formula 


k 



T 


where T and H are the matrices 


T 

H 



ds, and 
P dx dy. 


In these formulas, the following notation Is used: 


• P is the matrix defining the stress components a ■ { c^il » 

a22» 012^*^ terms of the vector B of stress^ 

degree 8-of -freed om 

o “ P B (Order P * 3 x 5) 

# N Is the 3x3 matrix of material constants, and In the case 
"" of a viscous Incompressible fluid Is 


N 


1 


10 0 
0 1 0 
0 0 1 


• Q Is the area of the element and F is Its boundary 


• R Is the matrix defining the surface tractions S ■ {Sj, 
So}^ In terms of the stress parameters B: 


S - R B (Order R - 2 x 5) 
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This corresponds to 


the Cauchy 


stress principle. 


S 


1 




■ '$ > S 
2x1 



S 


no • nP8 “ 

^ ^ ^ ^ ^ ^ 


Ai 

R ■ n P n • matrix of components of unit 
~ ^ normal to r 


e L is the matrix defining the boundary displacements v « 

~ {vi(s), V2 (s)}T in terms of the 8x1 vector q of nodal 

displacement degrees of freedom 

V • L q (Order L « 2 x 8) 


9.2 Q MATRIX CALCULATION 


We wish to compute an additional matrix Q that is similar to T except 
that it is of order 1x8: 


Q 

1x8 


f I"!’ "2' I 

n 1x2 2x8 


The computation of Q proceeds as follows: 
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where 


ORa;‘:':al 
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T~T 


- y <*3 - y's 


II 

n ’ 


*2 ~ *3 


ds - dy 


y y y 

u - (1 - — ) + 7- q . ^ " y <13 

^3 ^3 2 ^3 


/ ■ [<1 - “l * ^ <’■2 - ■‘S) 4 "S 


dy 


yj y, (X2 - X3) 

-q3 + -<l2 + 2 I3 


Along boundary 3-4 


f " f ("i « + “2 


v) ds 


where 


a - y < X 3 - + (y3 - Y2)‘ 


^4 " ^3 


n. 


*3 " *4 


(¥A)1 


'X - X- 


^4 ■^(x 3 - xj ^2 
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[‘-(P:)l»-(Pr) 


ds ■ ( ) dx 

*4 " *3 




*3 ” *4 


|1(.■(P;))..•(P;)..1] 


I dx 
(x^-XaT 


(y 4 “ yj) (y^ - yj) (x^ - x^) (x^ - x^) 


2 **4 2 *^2 2 ‘*5 


<lo - 


- 


Along boundary 4-1 


/ ■/ 


(n^ u + v) ds 


where 




1 £ ’ 2 £ 


‘’4' 




qr, ds 


- dy 
^4 
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- y <- y q* + ^ y qj) 

^4 



2 


‘’A 






4/ M3 73M4 ■ -3M5J 


9.3 THE NEW PENALIZED STIFFNESS MATRIX 


We are no« ready to compute the new stiffness matrix for the Incom' 
presslble visc<Mis flow problem. 


Step 1: Compute the usual stiffness matrix K using N 



I 

3x3 


K' 


T -1 
T H T 


Step 2; Let A^ denote the area of the element. Using the Q matrix 
discussed earlier, compute the perturbation stiffness matrix 


- -7 • Q 

^ e 


Step 3: Add K and K to get the element stiffness matrix 

^ *w p 


K® + K 
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The theory behind these cslculatlons is given In Appendix A. 


These modifications have been Incorporated into the SPAR code, 
pages to the SPAR user's guide are presented In Appendix B. 
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10. RESULTS 

Solutions were attempted for four two-dimensional example problems. 

The first three consisted of viscous flow of an incompressible fluid within 
rigid boundaries. The fourth problem consisted of two elastic plates 
coupled by a Stokesian fluid. Reasonable results were obtained for the 
velocity field for the first three problems. Descriptions of the example 
problems with sketches, finite element grids, SPAR input data listings, and 
tabulated results, are presented in Appendixes € through E. A pressure cal- 
culation r<Hitine was implemented but calculated pressures were not reason- 
able and no results are presented in this report. 

The first example (Appendix C) consists of parallel flow through a 
straight channel with uniform pressure boundary conditions applied at the 
entrance and exit of the channel. This problem has a linear pressure dis- 
trlhuLl.in in the flow direction and a parabolic velocity profile in the 
transverse direction. An 8 x 8 element mesh was used. A plot of the finite 
element grid, with the transverse scale enlarged for clarity, is shown. As 
can be seen from the tabulated results, the velocity remained essentially 
constant in the flow direction verifying that incompressibility is enforced. 

The second example (Appendix D) is a plane slider bearing lubrication 
problem. This problem consists of a moving guide surface separated from a 
stationary slide block by an incompressible viscous lubricant. As can be 
seen from the dimensions on the sketch, this model has a length-to-wldth 
ratio of 900. This caused a problem with element aspect ratios In SPAT? with 
an 8 X 8 mesh, so an 8 x 18 mesh was used for this problem. Again, a plot 
of the model with the y-dlrecL.Vn scale enlarged for clarity is shown. The 
computed '"elocitles agreed very well with the analytical values. 
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The third exevple (Appendix E) conslete of Incompreeelble flow In a 
driven cavity. The problen consiata of a aquare box encloaed on three aldea 
containing a vlacoua Incoapreaalble fluid driven on the upper aurface with a 
uniform velocity. 
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OF FOOrt 


Appendix A 


The following equations govern the steady uniform flow of a viscous 
incompressible fluid: 


o'- • “ p,. ■ f* 

iJ.J *^‘1 1 


d. . = k (v. . + V. .) 

ij i»J J.i 


‘‘ij ■ TT°ij 


dj^j^ - div V ■ 0 


n, (a.". - p 6 - .) ■ S. on P 
J ij iJ I a 


V . 
L 


V. on r 
X V 


1 




in n 


(A.l) 


Here are the deviatoric stress components. 


ij ij 3 “ij kk 

p is the hydrostatic pressure, fj^ the components of body force per unit 
volume, the components of the deformation rate tensor, the 

velocity components, pi the viscoscity of the fluid, and and 

denote prescribed tractions and velocities on portions and of 

3 

the boundary F of the domain C R “ich unit exterior normal n . 
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OF POOR QuALi jY 


We shall monencarily sat f ■ 0 , s • 0 without loss in generality. 

We next introduce a special complenentary energy functional 4 defined 
on a space of self-equilibrating deviatoric stresses and hydrostatic pres- 
sures: 

Z - {(o', p) ^ ^ 1 j ■ P'i • 0* (A. 2 ) 


(o'. “ P 6* .) n. • 0 on F } 
“ij ij j a 


♦ ( o', p ) 




o', dx 
ij 



-P Sij) 


V 


0 . 
V. ds 
1 


(A. 3) 


Here 5 and P are spaces of stresses and pressures, respectively, 
defined on the closed body Q which contain functions sufficiently smooth 
chat the functional is well-defined. The functional $ is essentially 
that introduced by Bratineau and Atluri^. 

Formally, the Euler-Lagrange equations corresponding to the stationary 
condition 


V (a' , F) ; 6 (o', p) , (o' , o)) - 0 


are 

''(i,j) 


d. . 
ij 




V. . V V » 0 

11 


I 

J 


(A.4) 


(A. 5) 


Bratineau, C. , Ying, L. A., and Atluri, S. N., "Analysis of 
Stokes Flow by a Hybrid Method," Finite Element Flow Analysis , University 
of Tokyo Press, pp. 981-988, 1982. 
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guAUtV 


Our next step Is to Introduce a perturbation of ^ associated 
with the hydrostatic pressures p . Let e denote an arbitrary positive 
number. Then we define 



Then, for any fixed p^ & P > the functional o' 
differentiable, and coercive: 


(A.6) 




(A. 7) 


♦ (o'*,p) is concave, 


& ^ 




■ IT iis'ii" * * ‘=*'’0’ r' 

ITT Ill'll! * IMIo llAvll^ 


+ C(Pq, v) 


i.e. , 

♦g (o' , )-►-<» as 1 1 o'*! I -► + 00 

Likewise, for any fixed a'*C. 5 , the functional p (a^ p '/ is concave, 

differentiable and coercive ($ (o' , p )-*•*<*» as lip!! ®) . Hence, we 

e -o lilt 

can conclude that 

• for every e > 0, there exists a unique saddle point (o' , p^) 
of ♦ 

e 

• as e 0 , (ct' , P ) conveys to a critical point (o',p) of 

•^e c 

the functional ^ . 


A-3 


LOCKHEED-HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC TR D867285 


OF POOR QUALITY 


Notice chat (a' » P^) latisfies the variational aquations 


( “ -s — 0 ♦ V , . . V ) 0 . dx ■ 0 
2y ij (i»j) iJ 


i 

I ^"'^i.i ♦ ePg) p dx - 


(A.8) 


for all (a . p) ^ Z 


Thus, 

Pg • ■ ^ (A. 9) 

It appears that the use of is equivalent to appending Che complimentary 

energy with an exterior penalty term corresponding to the incompressibility 
constraint. 


We can relax the constraint a'. . 'P>- 

iJ.J 1 


0 (equivalently, (a',p) 


€ 5 X P) by introducing the functional. 


L :Ax5xP'*R 

c 


L (X , (o', p)) - (o', p)- + r A. (oj*. . - p,£> dx 

e - - e ~ I ^ ij*j *■ 

Jo 


2 3 

with A - (L (n)) , the Euler equations of which are (formally). 




V. . 

i.i 


'P’i 


^i 


0 

0 

0 

0 


(A. 10) 


i.e., X is the velocity field defined on Che interior of the ■•cmain H . 
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OF POOR QUAliTY 


Lae ui now consider Che consCrucr i^^n of an astumed-scress hybrid 
finite element ep^roximetion using the funccionel (and ultimately 
L^ . We begin in the tradition?! way by introducing approximations of 
the stress a' end the boundary velocities v form 


11 


'22 


12 


> - p (x) e 


.1 


v,<;) 


b ■ L q 


(A.ll) 


trhere P is a matrix of polynomials in local element coordinates Xj, r. 2 » 

B is a vector of stress parameters, L is a matrix of polynomials in 
Xjj X 2 and q is a vector of nodal velocities associated with boundary 
nodes. Likewise, we approximate the element hydrostatic pressure field 
P by 

p “ A(x)p (A. 12) 


where A is, again, a new matrix of polynomials of ^ is a vector of 

m 

pressure degrees-of-freedom. At this stage, p and A should be selected 

T 

so that o'. . - P. . (or, in matrix notation, D o ~ 7 p ■ 0). 

ij,ji 

With these notations, the functional for a typical finite element 

becomes 


(B . P) 


IT T T 

-—rnS ♦ B Tq - p"Qq 


♦ p M P 

V mm m 


(A. 131 
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where 


OF POOR QUALITY 


H - 

T - 

Q - 

M - 


J_ 

2u 


r p"p 

~ - 

f P^ n L d$ 

Jjt MM 

r n L di 

Jr - - ~ 

f A^ A dx 
JQ - - 


dx 


The discrete functions! sssuaes s stationary value wb^enever 


(A. 14) 


3« 

- 0 

30 - 


3« 


e ■ 0 


Thus 


- H 8 ♦ T q 


- Q q e M 


h”' T q 


T?'' S3 


and, therefore, 


• T 3 3 


where K is the perturbed stiffness laatrix. 


T 1 T - 1 

K -T H T- — QM 0 
~ e--- 


(A.15) 


(A. 16) 
(A. 17^ 


(A. 18) 


(A. 19) 


The matrix H ' T is the usual assumed^stress hybrid stiffness matrix 
for the element whereas -Q M Q is a penalty-type matrix associated 
with the constraint div v ■ 0 . Notice that the hydrostatic pressure has 
been eliminated completely from the formulation, and is computed a posteriori 
by the formula 
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c- 


p 


■ A p 


— A Q q 
c * - - 


(A.20) 


Inttrior Velocity Foraulatiou . An tnrichcd «pproxiiucion, which nuty 
lt«d CO coordiuACt-invariant: tciffnata macricca, la obcainad if wc rapcac 
cho abova calculaciona uaing cha parcurbad Lagrangian of (A.IO). TF.18 

nacaaaicacaa chac va incroduca an indapandaoc approxinacion of Cha incarior 
valocicy X - V of cha cyna 


X 



(A.21) 


where b is a macrix of "bubble functions" (ganarally vanishing on F) 
associaced with dagraa-of- freedom paramatars p . For a typical element, 
wa have. 


» P' ® 6 “ c p 


wherein (8 . p) is given by ( 13} and 


B - r b^ P dx 
‘ ' 

c ■ r b^ 7^ A dx 
Jfi ~ - ~ 


Thus, instead of (IS) we have 




(A. 22) 


“HB+Tq+By - 0 


-Qq + eMp-c y “ 0 


B B - c p 




J 


(A. 23) 
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from which we compute 


OF POOR QUALiiY 


!o 3 


T ?i 3 


(A.24) 


where 


h"4 

-0 -e — ~ 


S| a r”‘ c Q 




-IT 1 -IT 
R — cM c 

-e e — 


(A. 25) 


Hence , 


e - H (T + B S) q 


P - (Q + c^ S) q 


(A. 26) 


with S > Sq - c ^ S| . Finally, the element stiffness matrix is 


- (T + B^ S)^ h“' (T ♦ B^ S) - e"' (Q f c^ S)^ 


h’* (Q + S) 


(A. 27) 


and the hydro* catic pressure is ^iven by 


P --A p 


(A. 28) 


where p is defined in (A. 26). 
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Appendix B 


Included as an attachment to this appendix are update pages to the SPAR 
Structural Analyals System Reference Manual (NASA CR 158970-1) dated 
December 1978. These updates describe the use of a new processor, EKSF, 
used to generate element Intrinsic stiffness matrices for Incompressible 
viscous flow analyses. Also described Is the velocity vector version of 
PLTB, PLTB/WEC, used for plotting flow vectors I..ulcatlng the magnitude and 
direction of velocities. 



Attachment to Appendix B 


Update pages to the SPAR Structural Analysis Syston 
Reference Manual (NASA CR 158970-1) 


Sect ion 


3.2 ELD- ELEMENT DEFINITION PROCESSOR 

3.2.1 General Rules, ELD Input 

3. 2. 1.1 Error Conditions 

3. 2. 1.2 Element Reference Frames 

3. 2. 1.3 Element Group/Index Designation 

3. 2. 1.4 The M)D Command 

3. 2. 1.5 The INC Command 

3.2.2 Structural Element Definition 

3. 2. 2.1 Line Elements 

3. 2. 2. 2 Area Elements 

3. 2. 2. 3 Three-Dimensional Elements 

3.2.3 Thermal Eleti«nt Definition 

3.2 E- E-STATE INITIATION 

3.4 EKS- ELEMENT INTRINSIC STIFFNESS AND STRESS MATRIX GENERATOR 

3.5 EKSF- INCOMPRESSIBLE VISCOUS FLOW ELEMENT INTRINSIC STIFFNESS 

MATRIX GENERATOR 

4 SPAR FORMAT SYSTEM MATRIX PROCESSORS 

4.1 TOPO- ELEMENT TOPOLOGY ANALYZER 

4.2 K- THE SYSTEM STIFFNESS MATRIX ASSEMBLER 

4.3 M- SYSTEM CONSISTENT MASS MATRIX ASSEMBLER 

4.4 KG- S'ibTEM INITIAL STRESS (GEOMETRIC) STIFFNESS MATRIX ASSEMBLER 

4.5 INV- SPAR FORMAT MATRIX DECOMPOSITION PROCESSOR 

4.6 PS- SPAR FORMAT MATRIX PRINTER 

5 UTILITY PROGRAMS 

5.1 AUS- ARITHMETIC UTILITY SYSTEM 

5.1.1 Miscellaneous 

5.1.2 General Arithmetic Operations 

5. 1.2.1 SUM 

5. 1.2. 2 PRODUCT 

5. 1.2. 3 UNION 

5. 1.2.4 XTY, XTYSYM, STYDIAG 

5. 1.2.5 NORM 

5. 1.2.6 RIGID 

5. 1.2. 7 RECIP, SORT, SQUARE 

5. 1.2. 8 RPROD, RTRAN, RINV 

5. 1.2. 9 LTOG, GTOL 

PREC6DJNG PAGE BLANK NOT RLMED 

iii 


Table 1-2: SPAR PROCESSOR FUNCTIONS 


NAME AMD 

SECTION 

REFERENCE 

TAB 3.1 


ELD 3.2 


E 3.3 


EKS 3.4 
EKSF 3.5 
TOPO 4.1 


K 4.2 

M 4.3 

KG 4.4 

FSM 12 


INV 4.5 


FUNCTION 

Translates user Inputs Into data sets containing basic 
tables of Information such as: 

- Joint Locations 

- Material Constants 

- Element Section Properties 

- Joint Reference Frame Orientations 

- Constraint Conditions 

- Rigid Lumped Mass Data 

(See Section 3.1 for a complex list) 

Produces data sets containing basic element definitions, 
l.e., connected lolnts. Integers pointing to applicable 
lines In tables of section properties, material constants, 
etc. 

Generates a system of data sets called the ’E-state, ' 
consisting of individual element Information packets 
containing data such as element geometry (dimensions, 
orientation), and literal section properties. E also forms 
the system diagonal mass matrix. 

Computes element stiffness and stress Influence matrices, 
and Inserts them Into the ’E-state.' 

Computes incompressible viscous flow "stiffness" matrices 
and inserts them Into the 'E-state.' 

Analyzes element Interconnection topology, and produces 
data sets used to guide other SPAR processors In forming 
and factoring assembled system matrices. 

Forms system elastic stiffness matrix. 

Forms system consistent mass matrix. 

Forms system geometric (pre-stress) stiffness matrix. 

Forms system matrices (dllltatlonal strain energy, 
gravitational energy, kinetic energy) associated with fluid 
elements. 

Factors system matrices In SPAR's standard sparse-matrlx 
format, e.g., K, K+KG, K-CM. 


1.2-4 


Section 3 

STRUCTURE DEFINITION 


To define the basic finite element model of the structure, the user 
proceeds as follows. 


- Execute TAB to define joint locations, joint reference frame 
orletitatlons, tables of section properties, and other basic 
components of the problem definition, as sumnarlzed on Table TAB-1 
In Section 3.1. 

- Execute AUS/TABLE to generate tables of section properties for 
three-dimensional solid and fluid elements, if required, as decrlbed 
In Section 3. 2. 2. 3. 

- Execute ELD to generate data sets containing basic element 
definitions, l.e., connected joints. Integers pointing to applicable 
lines In tables of section properties, etc. 

- Execute E to generate a system of data sets called the "E-state," 
consisting of Individual element Information packets containing data 
such as element geometry (dimensions, orientation), and literal 
section properties. 

- E also produces the sysLem diagonal mass matrix. 

- EKS Is executed tc compute individual element stiffness and stress 
recovery matrices, and Insert them Into the E-state. 

or 

- EKSF Is executed to compute individual element Incompressible 
viscous flow intrln*- Ic "stiffness" matrices, and insert them Into 
the E-state. 

All of the basic structural definition data sets produced as outlined above 
should be retained In Library 1. 


3-1 


Table 1-1: SPAR ELEMENT REPERTOIRE 


Name 

E21 

E22 

E23 

E24 

E25 


E31 

E32 

E33 

E41 


E42 

E43 

E44 


S4l 

S61 

S8l 


F41 

F61 

F81 


Description 

General straight beam elements such as 
such as channels, wlde-f langes, angles, 
tubes, zees, etc. 

Beams for which the intrinsic 
stiffness matrix Is given 

Bar - Axial Stiffness only 

Plane Beam 


See Volume 
3.1.7 - 9 

3.1.10 

3.1.11 

3.1.12 


1 Sections 


Zero-Length Element Used to Elastically 3.1.10 

Connect Geometrically Coincident Joints 

Two-Dimensional (area) Elements 3.1.13 

Triangular Membrane 
Triangular Plate 

Triangular Combined Membrane and Bending 
El«sent 

Quadrilateral Membrane, or 2-D Incom- 
pressible Viscous Flow Element (when 
used with EKSF). 

Quadrilateral Plate 

Quadrilateral Combined Membrane and Bending 
Element 

Quadrilateral Shear Panel 3.1.14 

Three-Dimensional Solids 3. 2. 2.3 


Tetrahedron (Pyramid) 

Pentahedron (Wedge) 

Hexahedron (Brick) 

Compressible Fluid Elements; 12., 3. 2. 2. 3 


Tetrahedron (Pyramid) 
Pentahedron (Wedge) 
Hexahedron (Brick) 


Notes : 

See Section 7.2 for examples of stress output 

See Volume 2 (theory) for element formulation details 

Aeolotroplc constitutive relations permitted, all area elements 

Laminated cross sections permitted for E33, E43 

Membrane/bending coupling permitted for E33, E43 

E41, E42, E43, E44 may be warped 

Aeolotroplc constitutive relations permitted for 3-D solids 
Non-structural mass permitted for line and area elements. 



EKSF 




3.5 EKSF-INCOMPRESSIBLE VISCOUS I'LOW ELEMENT INTRINSIC STIFFNESS MATRIX 

GENERATOR 

Function. EKSF functions similarly to EKS, !.<»., based on the 
dlnenslnns, section properties, etc., currently embedded In the element 
Information packets originated by processor E, EKSF computes Intrinsic 
stiffness and stress matrices for all elements other than E41 elements 
(e.g., E21 elements) and inserts them Into the packets. For E41 elements, 
EKSF computes Incompressible viscous flow "stiffness" matrices and Inserts 
them Into the packets. 

RESET Controls. Two additional reset controls have been Incorporated 
into EKSF which apply to E41 elements only. 


RESET Controls 

Default 


Name 

Value 

Meaning 

ELIB 

1 

Library containing the element Information packets. 

TIME 

0 

Nonzero value causes printout of intermediate CP and 
wall clock times. 

GAZERO 

10. "20 

Zero-test parameter, (beam area) x (shear modulus). 

CIZERO 

10.-20 

Zero-test parameter, beam non-uniform torsion constant. 

EPSILN 

.001 

The penalty parameter, e, used to enforce 
Incompressibility. 

XMU 

10.-^ 

Shear viscosity of the fluid. 

Note: 

EPSILN and 

XMU apply to E4l elements only. 


Core Requirements. EKS requires only a buffer area through which 
element information packets are transmitted. About 5,000 - 15,000 locations 
are usually suitable. 10 counts will vary In Inverse proportion to core 
space. 


3 . 5-1 


PLTB 




Velocity Vector Version (PLTB/WEC) 

PLTB/WEC functions like PLTB except that for all "detormed" plots, the 
"displacements" (flow velocities when executing in the viscous flow mode) 
are plotted with arrows indicating the magnitude and direction of the joint 
"displacements." The control statement DNORM remains In effect for 
normalizing joint "displacements" (velocities). Options 2A and 25 do not 
apply to this processor. All "deformations" are plotted as flow vectors. 
(See examples in Appendix.) 

Note: PLTB/WEC is available for plotting on the FR-80 plotter only. 
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OKIGl-Ai. rACl 
OF f'OOR QUAU.T 


Table C-1 

INPUT DATA FOR MODEL PROBLEi-; A 


PEARjNBIN202*SP AR< 1 ) 

.TEST A2/R 




1 

(dXQT 

TAB 





2 

START 

81 

3 4 5 6 

% 

:-o , 

NO nOTATIONS 

3 

TITLE 

• MODEL PROBLEM 

A 



4 

>*ATC 






5 

1 

30 

.♦06 .33 




6 

JLOC 






7 

1 

C. 

0. 0. 

T. 

.01 

C. 9 : 

8 

9 

.08 

0. 0. 

. 08 

.01 

C. 

9 

SA 






10 

1 

1 . 

J 




11 

CON = 

1 





12 

ZERO 

1 

,2: 9 ,61 ,9 





13 ZERO 2: 1,73,9: 2,8: 74,60 

14 dXjT ELD 

15 E4 1 

16 1 11 11 2 1 8 6 

17 SXQT E 

18 .iXQT EKSF 

19 RESET XMU=2.-04 

2C *PMD,E 

21 mXaT TOPO 

22 «XUT K 

23 RESET SPCP ? 

24 iXCT INV 

25 axuT AUS 

26 StSVEC 

27 applied forces 

2 8 C-'SE 1 

29 HI: J = 2,b: 1.5 

:C I-l: J=l: .75 

31 1=1: J=74,e0: 1.02 

32 1=1: J=73: 0.51 

33 OXC.T 5S0L 

34 dXLT VPRT 

35 PRINT STAT REAC 1 J 

36 PRINT SI AT OISP 1 1 

37 iX&T DCU 

38 TOC 1 
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OK POOR QUALIKY 


Table C-2 

COMPUTED VELOCITIES FOR MODEL PROBLEM A 


STATIC DISPLACEMENTS. 


NT 

1 

2 

^b 

. 73 5*03 

-. 6b4 ♦nu 

1 

. 31 t^O** 

.HOC * 

27 

.000 * 

. 0 d 0 

7 

. 311<D4 

.TCO * 

28 

. 31 5^J4 

. non 

3 

.296^C‘* 

.000 • 

29 

. 31 L*Q4 

' 1 1 7^11 

<4 

.271^o4 

. 000 ♦ 

30 

. 295^U4 

. 192^P1 


. 236^04 

. ODD » 

31 

.270^04 

.225>ni 

6 

. 1914 j4 

. OCO ♦ 

52 

.23b4^4 

. 242^01 

7 

. 137>u4 

.non • 

3 3 

. 192^J4 

. 199^ni 

8 

. 732^d3 

.000 * 

34 

. 1 3 3 ♦ 4 

. 137*01 

9 

.con * 

.000 * 

35 

- 738> J3 

.bLC*Ou 

in 

. 315^J4 

.000 ♦ 

3b 

.OOC * 

.con 

1 1 

. 31 j4 

. 23n^ni 

37 

. 314^04 

. 000 

12 

.295v^4 

. 394^-^l 

38 

. 309>D4 

. 29 4 »n.T 

1 3 

.270^04 

. 46) ♦ni 

39 

. 295*G4 

.4a3«02 

1*4 

.236^04 

. 4m«4 ♦Ol 

40 

.270^04 

. 353*00 

15 

. 192*w4 

. 3*>7^01 

4 1 

.23b^Q4 

-. 3C5*00 

lb 

. 1 38>u4 

.21C-401 

•42 

. 1 42^04 

-. 578*no 

1 7 

. 738*o3 

. 7io^no 

43 

. ) ’e>04 

-.664*00 

18 

.oco * 

. "CO * 

44 

. 739*03 

43b*00 

19 

. 31 5^04 

.000 • 

4 5 

.OCC * 

.OCO 

.r 

. 31C >04 

~. 1 V2^C) 

4 b 

. 31 5>04 

n "1 r> 

r l u 

2! 

.2*’5* J4 

318 ♦Oi 

47 

. 31C^C4 

9fl6* 

22 

. 2 7 L ♦ J 4 

-. 3c3*r*l 

48 

. 295^04 

161*01 

2 3 

.23b> v4 

-. 353^P1 

49 

. 270^04 

-. 158*0i 

2*4 

. 1 92^04 

-. 277^01 

50 

.23b»34 

730*00 

25 

. 1 38^C4 

175^01 





(Continued) 
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OF FOOH 

Table C-2 (Concluded) 


£.! 

. 192404 

-.112400 

b2 

• 1 3 8 ♦ j 4 

. 337400 

53 

. 7384J3 

. 336400 

54 

.•300 ♦ 

. 3 u 0 ♦ 

55 

. 31 5454 

.000 * 

36 

.310404 

. 27<=40„ 

57 

.295404 

. 376400 

58 

’704J4 

. 14 7400 

59 

. t ^ 6 4vj4 

-.514400 

6? 

.192404 

-.505400 

ol 

. 1384J4 

-.783400 

o2 

. 737403 

-.434400 

63 

. COO ♦ 

• 300 * 

64 

. 31 54 34 

. 3QC « 

o5 

. 310404 

.467400 

66 

.295434 

.895400 

67 

.271404 

.1274?! 

68 

.236404 

. Ic24ni 

69 

. 192404 

. 153401 

70 

. 1 3 84 j4 

. 109401 

71 

. 737403 

.444400 

72 

.C"0 ♦ 

.002 4 

73 

. 314404 

.''0 0 * 

74 

. 309404 

. 000 4 

75 

.295404 

.330 4 

76 

.270404 

.303 ^ 

77 

.236404 

. 300 4 

76 

. 1924 j4 

. 330 4 

7*3 

. 136404 

. 0 3 w 4 

6C 

. 741403 

.nco 4 

S 1 

.OOC * 

. 303 4 

EXIT 

2.792 

0 
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PLANE SLIDER BEARING 
MODEL PROBLEM B 
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r:\:G 

OF POOR 


Table D-1 

INPUT DATA FOR MODEL PROBLEM B 


PEAHSNBIN202*SPAR ( 1 ).TESTB/R 

1 

oiXQT TAB 

2 

START 171 3456 

3 

TITLE* MODEL PROBLEM 

4 

MATC 

5 

1 30. ♦OB .33 

6 

JLOC 

7 

1 0. C. 0. 

8 

9 .36 G. Oa 

9 

SA 

10 

1 1.0 

11 

CONrl 

12 

ZERO 1,2: 9,171,9 

13 

ZERO 2: 1,163,9 

14 

NONZERO 1: 1,163,9 

15 

BXQT ELD 

16 

E41 

17 

1 10 11 2 1 18 

10 

uiXQT E 

19 

8XQT EKSF 

20 

RESET XMU=2 -u4 

21 

SXUT TOPO 

22 

8X0T K 

23 

RESET SPDP 2 

24 

8XQT INV 

25 

8X0T AUS 

26 

SYSVEC 

27 

APPLIED FORCES 

28 

CASE 1 

29 

1=1: J=2,3:1.G 

30 

1=1: J=164,170: 0-5 

31 

SYSVEC 

32 

APPLIED MOTIO*'. 

33 

CASE 1 

34 

1=1: J=l, 163,9; 100.0 

35 

5XQT SSOL 

36 

aXOT VPRT 

37 

PRINT STAT REAC 1 1 

38 

PRINT STAT OISP 1 1 

39 

aXOT DCU 

40 

TOC 1 


2-0 , NO rotations 


• 00 Ob o • 9 I 

•0004 Q. 


19 
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OF POOR QiJ^LiTV 


Table D-2 

COMPUTED VELOCITIES FOR MODEL PROBLEM B 
STATIC DISPLACEMENTS. 


JOINT 

1 

2 

1 

. 100*03* 

.000 

2 

. 781*02 

-.179*00 

3 

.622*02 

-.355*00 

•4 

.460*02 

-.413*00 

5 

.332*02 

-.464*00 

6 

.225*02 

-.401*00 

7 

.125*02 

-.335*00 

8 

.726*01 

-.163*00 

9 

.000 ♦ 

.000 

ID 

. 100*03* 

.000 

1 

.828*J2 

. 156*00 

12 

, 636*02 

. 3C7*00 

13 

. 4n6*02 

.354*00 

14 

. 346*02 

.395*00 

15 

.223*u2 

. 339*00 

16 

.133*22 

.283*00 

17 

. 363*ol 

.136*00 

18 

.000 ♦ 

.000 

19 

. 100*03* 

. ODD 

:o 

. 793*02 

-.139*00 

21 

. 643*02 

-.271*00 

22 

. 486*j2 

-. 313*00 

23 

. 36J*02 

-.350*00 

24 

.251*02 

-.302*00 

2F 

. 146*02 

-.252*00 


26 

. 846*01 

-. 123*00 

27 

.DDO * 

.000 

28 

.100*03* 

.000 

29 

.840*02 

. 117*''0 

30 

. 657*02 

.225*^0 

31 

. 514*02 

.255*00 

32 

. 374*02 

.263*00 

33 

.250*02 

.241*00 

34 

.155*02 

.200*00 

35 

.491* j1 

.953-01 

36 

.000 * 

.000 

37 

. 100*03* 

.000 

38 

.806*02 

-.101*00 

39 

.666*D2 

-. 195*30 

40 

.514*02 

-.222*30 

41 

. 389*u2 

-.246*00 

42 

.278*02 

-.212*00 

43 

.167*02 

-. 176*00 

44 

.971*01 

-.861-01 

45 

.000 * 

.000 

46 

.100*03* 

.000 

47 

. 853*02 

.803-01 

48 

.679*02 

.151*00 

49 

.543*02 

. 166*00 

50 

.405*02 

.181*00 


(Continued) 
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Table D-2 (Continued) 


51 

.279*02 

.151*00 

52 

.178*02 

.124*00 

53 

.629*01 

.574-01 

54 

,000 ♦ 

.000 

55 

.100*03* 

.000 

56 

.821*02 

-.673-01 

57 

. 691*02 

-.126*00 

58 

.545*02 

-.140*00 

59 

.421*02 

-.153*00 

60 

,308*02 

-.130*00 

61 

.191*02 

.107*00 

62 

.111*02 

-.521-01 

63 

.000 * 

.000 

64 

, 100*03* 

.000 

55 

. 8 6 8 *C 2 

.470-01 

o6 

.704*02 

.834-01 

67 

.575*02 

.365-01 

6 8 

.439*02 

.391-01 

69 

. 311*02 

.697-01 

70 

.204*02 

.538-01 

71 

.781* j1 

.221-01 

72 

.HOC « 

* 0 >J 0 

73 

.100*03* 

. 000 

74 

,837*02 

-.360-01 

75 

.718*02 

-.630-01 


76 

,579*02 

77 

.458*62 

78 

. 342*02 

79 

.218*02 

SO 

.127*02 

81 

• COO * 

82 

.100*03* 

33 

.885*u2 

84 

,733*J2 

85 

.611*02 

86 

.478*02 

87 

• 348*02 

38 

.234*02 

89 

.953*01 

90 

.000 * 

91 

. 100*03* 

92 

.856*02 

93 

.750*02 

94 

.618*02 

95 

. 499*Li2 

96 

. 381*02 

97 

.250*02 

98 

. 145*02 

99 

.000 * 

IJO 

♦ 10G*U3* 


-.671-01 
-.693-01 
-.573-01 
-.<♦ * 44-01 
-. 210-01 
.000 ♦ 
.000 ♦ 
. 165-01 
.22C-01 
.139-01 
.512-02 
-.501-02 
-. 111-01 
-.109-01 
.000 * 
.000 * 
-.719-02 
-.584-02 
-.757-03 
.670-02 
. 924-"2 
.130-01 
. 731-02 
.000 ♦ 
.000 ♦ 


(Continued) 
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lai 

1:3 

1Q4 

1Q5 

136 

107 

108 

109 

110 
111 
112 

113 

114 

115 
.^16 

117 

118 

119 

120 
121 
122 

123 

124 

125 
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Table D-2 (Continued) 


. 905^02 
. 768^Q2 
. 654^D2 
. 524 + r ’ 

• 391 ♦ij2 
• 268*32 
• 115*32 
■ 000 * 

• 100+03* 
. 878*02 
. 787*02 
.665+Q2 
. 549*j2 
.428* J2 
• 281*32 
.167*02 
.000 * 
.100*034 
.929*02 
. 809*02 
• 7D6+L2 
.580*02 
. 443* '2 
. 309*02 
. 139+J2 


-.115-01 
-. 344“'31 
-.529-01 
-.727-''! 
-.747-01 
-.719-01 
-.420-01 
.000 4 

.000 4 

. 193-71 
.466-01 
.600-01 
. 764-01 
. 7u2-01 
.655-01 
. 331-'"! 

. 000 4 

.000 4 

-.373-01 
-.667-01 

— • 1 1 s * n 0 

— . 1 t 6 * OU 
-.141*00 
-.129*00 
-.715-01 


126 

. C 0 U 4 

.000 4 

127 

.100*034 

.^00 4 

128 

.905*02 

.438-01 

129 

.834*02 

.950-01 

132 

• 723*02 

.116*00 

131 

.611*02 

. 141*00 

a32 

.486*02 

. 126*00 

133 

. 334*02 

.113*00 

134 

.194*02 

.562-01 

135 

.000 4 

.000 4 

136 

. 100*034 

.000 4 

137 

.959*02 

-.610-01 

138 

. 8r.2*u2 

-.135*00 

139 

. 7‘)2*C2 

-. 174*00 

140 

.649*02 

-.215*00 

iHl 

.508*02 

-.203*00 

142 

.361*02 

-. 164*00 

143 

.170*02 

-.998-01 

144 

.000 4 

.000 4 

145 

. 100*034 

.000 4 

146 

.939*02 

.664-01 

147 

.393*^2 

. 139*00 

148 

.797*02 

. 167*30 

149 

.691*u2 

. 199*00 

130 

. 561*02 

.176*30 


(Cent inued) 
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Table D-2 (Concluded) 


ISl 

• 394*02 

. l56-*00 

152 

,229* J2 

.765-01 

153 

• □no « 

.000 

15R 

• lDO+03* 

.con 

155 

. 100403 

-.827-01 

156 

.930402 

-. 180*00 

157 

.857*02 

-.229*30 

158 

.741*32 

-.260*00 

15^ 

• 594*02 

-.263*00 

loO 

. 430*02 

-.237*00 

161 

.209*02 

127*00 

162 

.000 ♦ 

.000 

163 

. 100*C3% 

. 000 

164 

. 986*J2 

.868-01 

165 

.973*02 

.179*30 

166 

.897*02 

.213*30 

’67 

.797*02 

.250*00 

168 

.660*02 

. 220*00 

169 

.473*02 

. 193*00 

170 

.276*02 

.932-01 

171 

.000 * 

.000 

IT 

9.477 

u 
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MODEL PROBLEM C 
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Fig. E~2 - Driven Cavity Problem 

Computed Velocity Field 
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cm 

OF 


POOR O'iAUTY 


Table E-1 

INPUT DATA FOR MODEL PROBLEM C 


P£ARbN0lN202*SPARn ), TEST C/P 

1 dXUT tap 

2 START 121 3 *♦ 5 6 S 


3 

TITLE* MODEL PROBLEM 

C 

H 

MATC 


5 

1 30. ♦OS .33 


6 

JLOC 


7 

1 0. G. Q. 

1 

8 

11 c. i.n 0. 

1 

9 

SA 


10 

1 1.3 


11 

CONrl 


12 

ZERO 1,2: 1,111,11 

13 

ZERO 1,2: 11,121, 

11 

14 

ZERO 1,2: 2,10 


15 

ZERO 2; 112,120 


16 

NONZERO 1: 112,120 


17 

aXQT ELO 


18 

E41 


19 

1 2 13 12 


20 

dXOT E 


21 

OXOT EKSF 


22 

RESET XHU=.l 


23 

dPMD,E 


24 

aiXOT TOPO 


25 

iXQT K 


26 

RESET SPDP 2 


27 

^XWT INV 


28 

tfXOT AUS 


29 

SYSVEC 


30 

APPLIED MOTIONS 


31 

CASE 1 


32 

1=1: J=112,12L: ICC 

.0 

33 

iXQT SSOL 


34 

jXCT VPRT 


35 

PRINT STAT REAC 1 1 


36 

PRINT STAT OISP I 1 


37 

oiXOT DCU 


38 

TOC 1 



C-0 , NO POTATIONS 


l.C D. 


11 1 


in 10 


11 
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Table E-2 

COMPUTED VELOCITIES FOR MODEL PROBLEM C 


STATIC OISPLACCMENTS. 


JOINT 

1 

2 




1 

.000 * 

.000 * 

26 

-. 524*01 

.332*01 

2 

.000 * 

.000 ♦ 

27 

-. 123*u2 

. 136*31 

3 

.000 * 

.000 * 

28 

-. 344*j1 

-.239-r,l 


.300 * 

. 000 * 

29 

-.122*02 

-.136*31 

5 

.ooc * 

.000 ♦ 

30 

-.512*01 


6 

.000 • 

.000 ♦ 

31 

-.657*01 

-.242*01 

7 

.000 * 

.000 ♦ 

32 

.588*00 

-.297*31 

8 

.000 ♦ 

.000 4 

33 

.000 4 

.DCC 

9 

.300 « 

.000 4 

34 

.300 4 

.030 

10 

.000 * 

• OJO 4 

35 

-.447*31 

.496*31 

11 

.000 * 

.000 4 

36 

-. 320*01 

.943*31 

12 

.000 ♦ 

.ooc 4 

37 

-.156*02 

.641*01 

13 

-.179^01 

-.215*10 

38 

-.127* : 2 

.426*01 

14 

-.951*U0 

.138*01 

39 

-.215*^2 

-.972-01 

15 

-.4 53*01 

. ?C2*0C 

4C 

-.126*02 

-.439*01 

16 

-.330*01 

.568*0C 

41 

-.156*32 

-.639*01 

17 

-.584*21 

-.337-01 

42 

309*01 

-.834*01 

18 

-. 326*Jl 

-.546*00 

43 

-.446*01 

-.4o4*0l 

19 

-.440*01 

-.235*00 

44 

.000 4 

.3CC 

20 

-.901*00 

-.134*01 

45 

.300 4 

.OCO 

21 

-.1 78* J1 

.218*00 

46 

.113*31 

.103*02 

22 

.DOC • 

. DUO 4 

47 

-. 129*32 

.138*32 

23 

.COO * 

.000 4 

48 

-.114*02 

. 140*02 

24 

.539*00 

.303*31 

49 

-.275*02 

.779*01 


-. 662*C1 

.246*01 

50 

-.202*02 

-.107*00 


(Cont inued) 
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Table E-2 (Continued) 


51 

-.275^a2 

-.795+01 

52 

113-H]2 

-.141+02 

53 

-.128*u2 

-.137+02 

54 

.1174Q1 

-.101+02 

55 

.COO ♦ 

.000 

56 

♦ 000 * 

.000 

57 

-.739>01 

. 146+02 

58 

444+01 

.226+02 

59 

-.247+02 

.220+02 

6C 

-.208+02 

. 140+02 

61 

- . 363+02 

-.722-01 

o2 

-.2C8+02 

-.142+32 

63 

-.247+02 

-.221+02 

o4 

-.450+01 

-.224+02 

t>5 

-.745+01 

-.144+32 

66 

.000 * 

.000 

67 

.000 ♦ 

.OCC 

68 

.243+01 

.215+02 

69 

-. 18S+U2 

.319+02 

7G 

-.129+02 

.291+32 

71 

-.3l9+;.2 

.201+02 

72 

-,205+02 

-.144+00 

73 

-.320+02 

-.222+^2 

74 

-.132+02 

-.290+02 

75 

-.189+02 

-. 319+02 


76 

.247 + 01 

-.214+02 

77 

.000 ♦ 

.000 

78 

♦ 3DG ♦ 

. 033 

79 

-. 125+32 

.296+^2 

6 0 

-. 351+01 

.380+02 

81 

-.269+02 

. 384+02 

82 

-. 719+Jl 

.120+02 

83 

-.127+^2 

.247-01 

84 

-. 749+Jl 

-.121+02 

35 

-.270+02 

-.384+^2 

66 

-. 361+01 

-.381+02 

87 

-.126+02 

-.296+02 

88 

.COO ♦ 

.000 

89 

.coc ♦ 

.030 

90 

.622+^1 

. 360 + ni; 

91 

-.283+32 

,532+02 

92 

. 786+01 

. 124+02 

9 3 

.126+02 

.115+02 

94 

.205+02 

.111+00 

95 

. 12 6 + 02 

-.114+02 

96 

. 776+31 

-.126+02 

97 

-.283+32 

-.533+02 

98 

. 630+01 

-. 379+02 

99 

.000 * 

. 0 u C ’ 

inn 

• W 

.000 « 

.000 


(Continued) 
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Table E-2 (Concluded) 


lu 1 

341*02 

.639*02 


1-2 

.287*02 

.999*00 


in3 

.435*02 

.117*02 


1-4 

.531*02 

-. 149*00 


105 

. 549402 

.123-'^! 


106 

.532*02 

.232*00 


107 

.436*02 

-.119*02 


i08 

.287*02 

-.107*01 


109 

-.342*02 

-.636*02 


ilO 

. 000 » 

.000 

♦ 

ill 

.000 * 

.000 

* 

^12 

.10Q*03<‘ 

• 0k> J 

* 

113 

. 100*03* 

.*100 

* 

114 

.100*03* 

.000 

* 

115 

.100*03* 

.000 

« 

116 

• 100*03* 

. 030 

* 

117 

.100*03* 

.000 

* 

lie 

.100*03* 

.000 

* 

119 

.100*03* 

.000 

* 

120 

.100*03* 

.000 

* 

121 

.000 * 

. 000 

♦ 


Exit 13.996 0 6 
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